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Abstract 

An estimate of the average freeze-out volume for multifragmentation events is pre- 
sented. Values of volumes are obtained by means of a simulation using the exper- 
imental charged product partitions measured by the 4tt multidetector INDRA for 
129 Xe central collisions on nai Sn at 32 AMeV incident energy. The input parameters 
of the simulation are tuned by means of the comparison between the experimental 
and simulated velocity (or energy) spectra of particles and fragments. 
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A better knowledge of multifragmentation properties is of the highest im- 
portance in the investigation of the liquid-gas phase transition in hot nu- 
clei [1,2,3,4]. In particular, in various statistical and thermodynamical ap- 
proaches, the concept of freeze-out volume is introduced, which can be de- 
fined as the volume occupied by the ejectiles of the multifragmenting source 
when their mutual nuclear interactions become negligible. Such a volume ap- 
pears as a key quantity [4] and its knowledge is particularly important in the 
extraction of fundamental observables such as the microcanonical heat ca- 
pacity and its negative branch or the shape of caloric curves under external 
constraints [1,5,6,7]. 

Up to now volume or density information at freeze-out was derived in var- 
ious ways. For example by comparing average static and kinetic properties 
of fragment distributions with statistical multifragmentation models in which 
the freeze-out volume is an input parameter [8,9,10,11] or from nuclear caloric 
curves using an expanding Fermi gas hypothesis to extract average nuclear 
densities [12]. In this work we present a more direct approach to determine 
freeze-out volumes. Indeed dynamical simulations show that a geometrical pic- 
ture is fully relevant on the event by event basis and can be used to estimate 
average volumes of a given class of events [13]. In that context we obtained 
values of the average freeze-out volume in multifragmentation events, from 
a "fused system" produced in central collisions, by employing a simulation 
directly built event by event from the data collected with INDRA [14]. At 
the present stage we do not want to have a fully consistent understanding 
of parameter values derived from simulations but rather a very good repro- 
duction of data using reasonable physical hypotheses. Further details of the 
experimental and calibration procedures may be found in Refs. [15,16]. 

Complete experimental events (total detected charge > 93% and total mea- 
sured momentum > 80% of the entrance channel values) with flow angle > 60° 
(corresponding to compact single source reactions [17,18,19]) for the reaction 
129 Xe+ na *Sn at 32 AMeV were selected. The requested completeness on the 
total detected charge, more severe than that usually employed by the INDRA 
collaboration (~ 80%) [15,16], was justified by the necessity of a freeze-out 
source as close as possible to the reality as input for the simulation, to avoid 
underestimations of the total Coulomb repulsion among fragments and parti- 
cles. Main properties of selected multifragmenting sources are summarized in 
figure 1. Depending on the required completeness on the total detected charge 
(> 93% and > 77%), average charged product multiplicity varies from 23.8 
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to 26.2 whereas fragment multiplicity (with charge > 5) increases from 4.13 
to 4.72. One can also verify that the largest completeness does not introduce 
substantial bias on relevant observables as the differential charge multiplicity 
distribution, the average experimental velocity of fragments or the width of 
their velocity spectra. Moreover flow angle distributions for both selections 
(not shown) exactly superimpose. 
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Fig. 1. Properties of selected collisions: (a) differential charge multiplicity distribu- 
tions, (b-c) total charged product and fragment multiplicities and (d) average and 
width (standard deviation in the inset) of the centre of mass velocity spectrum of 
fragments as a function of their charge, regrouped by two charge units. Open (full) 
symbols correspond to total detected charge greater than or equal to 77% (93%) of 
the entrance channel. Vertical bars are statistical errors. 

The multifragmenting source is reconstructed in the reaction centre of mass 
from all the fragments (Z > 5) and twice the particles [Z = 1,2) and light 
fragments [Z = 3,4) emitted in the range 60°-120° in order to partially ex- 
clude pre-equilibrium emission [20] . Fast isotropic emission can not be removed 
from the source without any theoretical assumption and, in that respect, cal- 
culated excitation energy, mass and freeze-out volume of the source should 
be considered as upper limits. Atomic mass of detected fragments (Z > 5) 
was calculated from the EAL formula [21]. The number of neutrons (which 
are undetected) was calculated to keep the N/Z ratio of the entrance chan- 
nel [15]. With such a procedure the average atomic mass and atomic number 



3 



Table 1 



Radii for light particles (H, He) and tq for light fragments, from [22] 



proton 


deuteron 


triton 3 He a 6 He 


Li 


Be 


B 


C 


1.03fm 


2.8fm 


2.2fm 2.4fm 2.2fm 2.4fm 


1.7fm 


1.5fm 


1.45fm 


1.3fm 



of the reconstructed sources are respectively 217 and 91. Then, the partition 
at freeze-out for each event was built by "dressing" detected fragments pro- 
portionally to their measured charge (calculated mass for neutron "dressing" ) 
with certain percentages of detected (calculated) particles, light fragments 
and neutrons. Those percentages constitute one parameter of the simulation. 
The dressed fragments, assumed to be spherical, and the remaining particles 
and light fragments, if any, were placed in a compact configuration with a 
minimum distance D min among the surfaces. It was realized by first putting 
particles and fragments at random on the surface of a big sphere with volume 
equal to 30 Vb (Vo being the volume of the source at normal density) and then 
moving fragments and particles one by one toward the centre of the sphere 
by homothetic steps. D min is another parameter of the simulation. The radius 
of each fragment was calculated according to the formula R = roAs, where 
A is the fragment mass and ro=1.2 fm; for light particles (light fragments) 
experimental radii (deduced ro) summarized in table 1 were taken. The radius 
for neutrons was chosen equal to the proton radius. 

For each event a consistent calorimetry was made to derive the excitation 
energy E* of the corresponding source which undergoes multifragmentation. 
Thus the partition between internal excitation energy (for fragments) and 
kinetic energy at freeze-out was determined. To do that a variation of the 
level density for fragments was introduced. The level density is expected to 
vanish at high excitation energies [23] and the formalism adopted here is that 
proposed in [24], where the level density at excitation energy e is expressed as 
the Fermi gas level density p FG modified by a modulation factor: 

p(e) = p FG (e)e~^ (1) 



This corresponds to introducing an intrinsic temperature for fragments Tf rag 
which verifies 



T frag 2<Kf°> T lim v ' 

where < K*° > is the average kinetic energy of fragments and particles at 
freeze-out and T Hm the maximum temperature attainable by fragments. T Um 
is a parameter of the simulation. 
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Equations used for calorimetry and to derive the sharing between internal 
excitation energy and kinetic energy on the event by event basis are the fol- 
lowing: 

M cp 

E* s +AB S = J2 K% + AB cp + Mf° < > +M^e frag + AB n (3) 
fc=i 

nfrag 

E* s + AB S = {Mf° - 1) < K'° > + £ a k 9% ag + AB fo + v£ ul (4) 

k=i 

where 0f rag is equivalent to the temperature Tf rag in an ensemble average. 
AB, s, K, cp, n, and fo stand respectively for mass defect, source, kinetic 
energy, charged products, neutrons and freeze-out. The neutrons evaporated 
from primary fragments, M^ vap , have an average kinetic energy along the de- 
excitation chain equal to the initial temperature of fragments 0f rag [25]. An 
internal excitation energy equal to Y^k=\ 9 a k^}ra g i s associated to the frag- 
ments at freeze-out, where a k = j^MeV' 1 [26], A k is the mass of the k th 
fragment and nfrag is the fragment multiplicity at freeze-out. is the to- 
tal multiplicity at freeze-out. The total kinetic energy (M^° — 1) < > is 
shared at random between all the particles and fragments at freeze-out un- 
der constraints of conservation laws (linear and angular momentum). Vq° u1 is 
the Coulomb energy of the space configuration for freeze-out previously deter- 
mined. Nuclear interactions between fragments or particles are neglected and 
we shall see later that this approximation is reasonable. A radial collective 
expansion energy Er can also be introduced in equations (3) and (4) accord- 
ing to the formula Er = J2tf=i (r^) 2 ^Eq, where Rq is the rms of fragment 
distances to centre at the freeze-out volume, E the radial expansion energy 
at R and r is the distance of the considered particle/fragment of mass A from 
the centre of the fragmented source. 

Particles and fragments were then propagated under the effect of their mutual 
Coulomb repulsion; during propagation fragments de-excited, by means of an 
algorithm largely inspired by the SIMON code [27,28]. The main differences 
concern the tuning of the emission time during the evaporation sequence in 
order to reproduce the results of theoretical calculations for neutron emis- 
sion [29] and the constraint that the evaporated particles are inside the list 
of the particles placed on each fragment at freeze-out. The used emission bar- 
riers come from experimental data ([30] for Z — 1,2 and [31] for Z = 3,4). 
In this way at the end of the de-excitation phase we obtain secondary charge 
(and mass) distributions for fragments close to the experimental (calculated) 
ones; final charges (masses) are recovered, within two charge units (four mass 
units), for 98% (85%) of fragments. Finally experimental angular and energy 
resolutions were taken into account. 
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The last step was the comparison between experimental and simulated spectra 
both for the energy of the particles and for the velocity of the fragments 
(average and width) as a function of their charge, in order to tune the four 
parameters of the simulation (E , D min , T Um and the percentages of evaporated 
particles). We chose to compare the velocity spectra instead of the energy 
ones because the velocity is less sensitive to the final mass of fragments at 
the end of the de-excitation process. The explored range for the percentage 
of evaporated particles was between 0% and 100% (no free particles at freeze- 
out); 30% was suggested by [15]. For D min we investigated an interval ranging 
from fm (maximum possible approach without overlap) up to 5 fm; for E we 
tested from (no collective radial energy) to 1.2 MeV. Finally the explored 
values for Tu m ranged among 6 MeV and 20 MeV; in previous studies on 
the same or similar sample of events [20,32], Tn m values in the range 10- 
12 MeV were derived. The limiting temperature influences mainly the width 
(standard deviation) of the velocity spectra; the percentage of evaporated 
particles controls all the studied observables (both the standard deviation 
and the average value of the fragment velocity spectra and also the energy 
spectra of light particles). The distance among the nuclei surfaces at freeze- 
out and the radial collective energy control mainly the average of the fragment 
velocity spectra and, more weakly, the particle energy spectra; D min and E 
are correlated, indeed a larger surface distance implies a weaker Coulomb 
repulsion which can be compensated for by a larger radial kinetic energy. 

A x 2 minimization procedure was used to determine the best fit to the data. 
To reduce the total influence of particles and light fragments but emphasize 
their high energy tails very sensitive to freeze-out emissions, x 2 was calculated 
using all the fragment velocity spectra and the Log of particle and light frag- 
ment energy spectra. The best agreement, corresponding to a x 2= 0.953, was 
obtained for the following set of parameters: T} irrt =10 MeV, 90% of evaporated 
particles and light fragments, D min =2 fm and E =0.Q MeV. The obtained val- 
ues for the percentage of evaporated particles are larger than those extracted 
in [15]; this discrepancy may be due to some intrinsic lack of sensitivity of the 
method proposed in [15]. Indeed fragment-particle correlations are only fully 
sensitive when fragments de-excite with a sufficient distance between them 
and in that sense it is impossible to go back up to the freeze-out configu- 
ration. In any case the values presented in [15] constitute a reliable inferior 
limit for the percentage of evaporated particles. The total average mass of 
dressed fragments is 208, which corresponds to 96% of the mass of the source 
at freeze-out. 

Note that with the retained parameters the calorimetry procedure gives an 
average excitation energy (thermal+collective) < E* >=6.7 AMeV for the 
source which undergoes multifragmentation, which leads, for fragments, to 
a temperature Tf rag =< 0f rag >=6.3MeV and an average excitation energy of 
3.9 AMeV. The average kinetic energy of fragments and particles at freeze-out 
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corresponds, if interpreted in terms of "kinetic temperature" , to a value T^n 
of 17.5 MeV. A deep understanding of those different numbers and their rela- 
tion with the observed fragment partitions is out of the scope of the present 
letter. However we can mention three possible explanations that one will have 
to discuss in a near future. One is fully related to the procedure followed here, 
assuming thermal equilibrium for fragments, which reveals the major role of 
a limiting temperature (for fragments); in that context the influence of Tu m 
on partitions in microcanonical multifragmentation model like [11] will have 
to be checked in details. A second one refers to a fast fragmentation for which 
particles and fragments are early emitted: average primary kinetic energies 
are then related to the Fermi momentum [33,34,35]. Finally the "kinetic tem- 
perature" deduced could also reflect the fact that few very energetic particles 
emitted during the expansion-thermalization phase [36] significantly increase 
the average kinetic energy related to T kin . 

In figure 2 the experimental centre of mass average velocity of the fragments 
(full symbols) is compared to the best simulation (open symbols) as a func- 
tion of the final fragment charge. From the inset it is possible to appreciate 
the small absolute gap between the experimental data and the simulation. 
Coulomb repulsion at freeze-out contributes to ~70-80% of the calculated av- 
erage velocities. The standard deviation of the fragment centre of mass velocity 
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Fig. 2. Experimental (full symbols), best simulated (open symbols) and bad simu- 
lated (line, see text) centre of mass average velocity of fragments as a function of the 
final fragment charge, regrouped by two charge units. Vertical bars are statistical 
errors. The inset shows the absolute difference between the experimental data and 
the best simulation. 
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spectra as a function of the fragment charge for the experimental data (full 
symbols) and the best simulation (open symbols) is presented in figure 3. Spa- 
tial configurations at freeze-out and fragment decays are only responsible for 
~60-70% of the observed widths and the introduction of a limiting tempera- 
ture in the simulation turned to be mandatory to account for the experimental 
values. In figure 4 the experimental and best simulated centre of mass energy 
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Fig. 3. Standard deviation for the centre of mass velocity spectra of fragments as a 
function of the fragment secondary charge, regrouped by two charge units (symbols 
and line as in figure 2). Error bars represent errors coming from the fits. 

spectra for particles and lithium fragments are presented. The agreement be- 
tween simulation and data is good. Relative angles and velocities between 
fragment pairs were also compared through relative velocity correlation func- 
tions (not shown) and exhibit a reasonable agreement between experimental 
and simulated data. 

Before obtaining estimates on freeze-out volume, we can discuss the error bars 
on the different parameters of the simulation. The standard method was used: 
errors are extracted from simulations for which \ 2 = bestx 2 + 1- The derived 
ranges for parameters are the following: 9-11 MeV for Tu m , 70-90% of evapo- 
rated particles and light fragments, 0.5-3 fm for D m i n and 0.3-0.6 MeV for Eq. 
As previously noted, only those two last parameters are strongly correlated. 
Therefore, the same reasonable % 2 is obtained for couples of extreme values 
like 0.5 fm-0.3 MeV and 2 fm-0.6 MeV keeping constant T Um and the percent- 
age of evaporation. Note that even in simulations with D min =0.5 fm average 
distances among the surfaces of nuclei equal to or larger than 2.4 fm are ob- 
tained, compatible with the freeze-out definition as the configuration where 
the nuclear force among the primary products vanishes. Finally to appreciate 
the sensitivity of parameters, as an example, results from a simulation with 
X 2 = bestx 2 + 3 are also displayed in figures 2, 3 and 4. For that simulation 
three parameters have values close to or at the limits of the error bars just 
discussed before, Tu m — 11 MeV, D min =0 fm and £'o=0.34 MeV, and one has 
a value largely outside: 0% evaporation. It corresponds to a compact freeze- 
out state containing all the detected charged products and calculated number 
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Fig. 4. Center of mass energy spectra per nucleon for light particles and lithium 
emitted with polar angle in the centre of mass 60° < i? cm < 120° (symbols as in 
figure 2 and histograms for bad \ 2 simulation). 

of neutrons (sphere volume of 7.4 Vo - see definition below) . We note a wrong 
slope in the spectrum of the average velocities of fragments versus their charge 
(figure 2) mainly related to the freeze-out topology imposed by the large num- 
ber of particles (the heavier the fragment, the larger its distance relative to the 
center of the volume). The width of the velocity spectra are underestimated 
(figure 3) due to the absence of smearing effect from evaporation. The energy 
spectra of particles and light fragments (figure 4) have steeper slopes than 
the experimental ones because of the absence of evaporation (no boost from 
primary fragments). 

Once the values of parameters and their range defined, the corresponding 
freeze-out volumes were calculated from the envelopes of nuclei at freeze-out 
(before starting the Coulomb propagation). The envelope was defined in such 
a way that two external adjoining spherical nuclei Ai and Aj (whose centres 
are located at distances di and dj from the centre of the fragmenting source) 
are linked by a portion of sphere with same centre but a radius R e = 1/2 
(di + dj + r Ay 3 + r Ay s ). It is worthwhile to recall that the volume is in- 
fluenced by two parameters of the simulation: the minimum surface distance 
D min and the percentage of evaporated particles, and by the compact con- 
figuration built at random. The technique used to calculate the volume of 
the envelope was to fill up uniformly with points a sphere including all the 
nuclei and to calculate the ratio among the number of points included inside 
the envelope and the total number of points inside the sphere. In this way it 
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is possible to obtain a good representation of the effective volume occupied 
by the nuclei also for non-spherical configurations. For the best simulation 
the average envelope volume is 4.2Vo; the volume of the sphere including all 
the nuclei is an overestimation and, on average, for the best simulation, it is 
7.65Vo. Volume estimates, taking into account error bars on parameters, are 
presented in figure 5. Envelope volumes range from 3.2 to 5.2Vo and sphere 
volumes from 5.7 to 9.6Vq. A small increase (~10%) of envelope volumes is 
observed when the fragment multiplicity increases from 2 to 8 whereas sphere 
volumes keep constant whatever the multiplicity. Standard deviations of en- 
velope volumes generated by the simulation realized event by event are close 
to 10% of average volumes, which qualitatively agree with predictions, in the 
coexistence region, of a microcanonical lattice gas model with a constrained 
average volume [7,37]. 




Envelope volume 



Sphere volume 
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V/V 

Fig. 5. Estimated average freeze-out volume in units of Vo- Envelope volumes refer 
to realistic volumes (see text) and sphere volumes correspond to spheres including 
all the nuclei at freeze-out. 

To summarize, with the powerful 4n multidetector INDRA, we were able to 
obtain information on the freeze-out volume in single source multifragmenta- 
tion events in the framework of a simulation using the experimental data. The 
presented method needs data with a very high degree of completeness, which 
is crucial for a good estimate of Coulomb energy. The use of the widths of frag- 
ment velocity spectra in the comparison between data and simulation shows 
that the introduction of a limiting temperature in the range 9-11 MeV seems 
mandatory. Work is in progress to study the evolution up to 50 AMeV incident 
energy of freeze-out volume for the same system and to derive consistent and 
reliable information from parameters of simulations [38]. 
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